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The temperature dependent stability of the magnetic phases of FeRh were investigated by means of 
total energy calculations with magnetic disorder treated within the uncompensated disordered local 
moment (uDLM) approach. In addition, Monte Carlo simulations based on the extended Heisenberg 
model have been performed, using exchange coupling parameters obtained from first principles. The 
crucial role and interplay of two factors in the metamagnetic transition in FeRh has been revealed, 
namely the dependence of the Fe-Fe exchange coupling parameters on the temperature-governed 
degree of magnetic disorder in the system and the stabilizing nature of the induced magnetic moment 
on Rh-sites. An important observation is the temperature dependence of these two competing 
factors. 

PACS numbers: Valid PACS appear here 


INTRODUCTION 


FeRh with composition close to being equiatomic crys¬ 
tallizes in the CsCl structure and exhibits rather inter¬ 
esting magnetic properties attractive for investigation for 
various reasons. It is antiferromagnetically (AFM) or¬ 
dered in the ground state and reveals a first-order meta¬ 
magnetic transition to the ferromagnetic (FM) state at 
Tm ~ 340 — 350 K [I|. A transition to a paramagnetic 
(PM) state occurs at Tc = 675 K [ij. It is worth noting 
that the temperature of the metamagnetic transition is 
very sensitive to the conditions of sample preparation. 
Despite several attempts to shed light on the physical 
origin of the magnetic prcmerties of FeRh-based alloys, 
both within experimental [2^ and theoretical [Ml in¬ 
vestigations, they are still under debate. 

Various mechanism have been suggested to explain the 
driving force for the AFM-FM transition. The early 
model suggested by Kittel 0, namely the exchange- 
interaction-inversion model, associated the AFM-FM 
phase transition with the dominant role of the change of 
magnetoelastic energy. More recent investigations, how¬ 
ever, demonstrated a minor role of the exchange magne- 
toelastic energy [l5|. 

The experimentally observed large change of the en¬ 
tropy, i.e. 14.0 mJ/g/K or 12.58 mJ/g/K [ll, is 
mainly attributed to the electronic contribution related 
to spin fluctuations on the Rh atoms. This observa¬ 
tion implies a key role of the Rh induced mag netic mo¬ 
ments for the stabilization of the FM state This 

idea was supported by theoretical investigations based 
on first-principles electronic structure calculations Si- 

0, E nil. 


So far there is no clear understanding of the finite 
temperature magnetic properties of FeRh when small 
amounts of impurities are present [HIIMII. Recent 
first-principles investigations by Staunton et al. [l3| . 


based on the analysis of the electronic entropy highlights 
the impact of small compositional changes on the tem¬ 
perature of the metamagnetic transition. 

The reversibility of this transition and the small relax¬ 
ation time makes it attractive for applications. Transport 
measurements demonstrate a strong drop of the electri¬ 
cal resistivity during the metamagnetic phase transition 
Til 2^ from the AFM to the FM state. As the meta¬ 


magnetic transition can be manipulated by an external 
magnetic field, this feature of the resistivity leads to a gi¬ 
ant magnetoresistance phenomena in the system around 
Tm, that makes FeRh an appealing material for future 
data storage devices Esin. 


COMPUTATIONAL DETAILS 

Within the present study spin-polarized electronic 
structure calculations have been performed using the 
fully relativistic multiple scattering KKR (Korringa- 
Kohn-Rostoker) Green function method lUj. All 
calculations have been performed in full-potential (FP) 
mode. Density functional theory employing the Gener¬ 
alized Gradient Approximation (GGA) was used with 
the parametrization of the exchange-correlation poten¬ 
tial as given by Perdew, Burke, Ernzerhof (PBE) [j^ . 
For the angular momentum expansion of the Green func¬ 
tion a cutoff of = 3 was applied. For determining 
configurational averages in substitutionally disordered al¬ 
loys, the self-consistent coherent potential approximation 
(GPA) method was employed. The magnetic disorder in 
the DFT calculations for systems in the paramagnetic 
(PM) state (T > Tc) was treated within the Disordered 
Local Moment (DLM) method. To simulate the tem¬ 
perature induced partial magnetic disorder below the 
Gurie temperature, the so-called uncompensated Disor¬ 
dered Local Moment (uDLM) approximation (see, e.g. 
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Ref. 27|) was used. In this case the Fe subsystem was 
represented by the pseudo-alloy Fe|_^Fe;^ with Fe sites 
occupied by two-components Fe^ and Fe"*" with the op¬ 
posite directions of magnetic moment, ’up’ and ’down’ 
respectively, and x varying in the interval [0.0,0.5]. 

The finite temperature magnetic properties have been 
investigated via Monte Carlo (MC) simulations based 
on the extended Heisenberg model, using a standard 


Metropolis algorithm [26 


The exchange coupling 


parameters for these calculations were obtained us¬ 
ing the expression given by Lichtenstein 0,|3ll. 


RESULTS 

The calculated total energies of the FM and AFM 
states of FeRh as a function of the lattice parameter, 
represented in Fig. [U are in full agreement with the ex¬ 
perimental results as well as with calculations of others 
authors The energy minimum for the AFM state 

for the lattice parameter a = 5.63 a.u. is about 17 meV 
lower than that for the FM state occurring at a = 5.66 
a.u. Thus, the AFM-FM transition should be accompa¬ 
nied by a lattice expansion with a magnetovolume effect 
of 1.6 % which slightly overestimates the value observed 
experimentally, 1 % [l^ . 



FIG. 1. Total energy calculated for the FM and AFM states 
of FeRh as a function of the lattice parameter a 

The Fe magnetic moment obtained for the energy min¬ 
imum in the FM state is 3.3 /rs, which is in line with the 
neutron scattering measurements giving 3.0 hb per Fe 
atom [^ . Its magnitude remains almost unchanged in 
the AFM ground state (3.2 /r^). The induced total mag¬ 
netic moment on the Rh atomic sites calculated for the 
FM state is equal to 1.0 fiB and vanishes in the AFM 
state. 

The density of states (DOS) curves for the FM, AFM 
and DLM states are shown in Fig. [2] (a), (b), (c), respec¬ 
tively. Only a weak dependency on the volume is observ¬ 
able for the DOS calculated for the FM (a) and AFM 
(b) states. In the FM state a strong spin-dependent hy¬ 


bridization of Rh states is apparent ([2] (a)), which leads 
to the formation of a magnetic moment of 1 /is on the 
Rh atomic site. In the AFM state the Rh DOS for the 
majority- and minority-spin states are identical due to 
symmetry, resulting in a Rh magnetic moment equal 
to 0 ^B- The Rh-related electronic states for both of 
the DOS spin channels in these cases exhibits hybridiza¬ 
tion with Fe minority-spin states (essentially, above the 
Fermi energy) and majority-spin states (below the Fermi 
energy), in line with the discussion by Sandratskii and 
Mavropoulos Q- This demonstrates that in the AFM 
state the vanishing total magnetic moment on the Rh 
site is a result of the hybridization-governed redistribu¬ 
tion of spin density and not because of vanishing spin 
density within the Rh atomic site. The DOS 

Note also, that the Rh DOS at the Fermi level is rather 
small and, therefore, a pronounced Stoner enhancement 
of the magnetic moment induced on Rh can be ruled 
out. Thus, as will be also discussed below, the large Rh 
magnetic moment in the FM state should be attributed 
to a strong spin-dependent hybridization. 

As the AFM-FM transition occurs at finite tempera¬ 
ture it is necessary to take into account the temperature 
induced magnetic disorder in the system when compar¬ 
ing the total energies of the two states. For this reason 
the calculations have been performed accounting for mag¬ 
netic disorder treated within the uDLM approximation. 
In the case of the FM state with partial magnetic dis¬ 
order, the normalized magnetic moment at each Fe site 
is mpe =< Alpe > /Mpe = (Mpe(l “ x) - Mp^x) /= 
(1 — 2x). In the case of the partially disordered AFM 
state, the same procedure was applied to each Fe sublat- 
tice, Mpg and Mp^ having opposite alignment of magnetic 
moments with respect to each other, i.e. fhp^ =< Mp^ > 
/Mpe =< Mpe, > /Mpe = (1 — 2x). Thus, in this case 
the total magnetic moment is equal to 0 for each rfip^ 
value. 

Fig. El shows the total energy difference Eapm — Epm 
as a function of rfip^ for FeRh with the lattice parameter 
a = 5.63 corresponding to the energy minimum of the 
AFM state. At a; = 0, the difference is negative, demon¬ 
strating the stability of the AFM state in line with the 
results shown in Fig. [TJ An increase of the disorder rep¬ 
resented by a decrease of (1 — 2a;) leads to a decrease of 
stability of the AFM state such that at rhp^ < 0.8 the 
FM state becomes more stable up to the fully disordered 
state with fhp^ = 1 — 2a; = 0., when both types of mag¬ 
netic order, FM and AFM, have the same energy. The Rh 
magnetic moment in the case of FM order exhibits almost 
a linear dependence on rfip^ changing from Mpp = 0.0 ^b 
in the fully disordered DLM state to = 1.0 ^jlb in 
the ordered FM state. 

In summary, these results demonstrate the following 
effects of increasing magnetic disorder: (i) a stabilization 
of the FM state with respect to the AFM state and (ii) 
the stability of the FM state is a result of the decrease of 
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FIG. 2. Density of states for the FM state (a), AFM state 

(b) and random-spin state treated within the DLM approach 

(c) . The solid and dashed lines in (a) and (b) represent the 
results for FeRh with the lattice parameters corresponding to 
the minimum of total energy for the AFM state (solid line, 
a = 5.63 a.u.) and AFM state (dashed line, a = 5.66 a.u.) 


the Rh magnetic moment. 

Therefore, when discussing the driving forces behind 
the metamagnetic AFM-FM phase transition, other ad¬ 
ditional effects than just the magnetization of the Rh 
sublattice have to be considered. This has alre ady be¬ 
come apparent within various investigations U M- 

Below we will investigate the features of the interatomic 
exchange interactions to demonstrate their crucial role 
for the AFM-FM phase transition. 

First we discuss some features of the Rh magnetic 



FIG. 3. Top: Difference between the total energies, Sapm — 
Efm as a function of rhpe, obtained by uDLM calculations 
for the FeRh compound having FM and AFM magnetic or¬ 
der, with lattice parameters corresponding to the total energy 
minimum in the AFM state, a = 5.63 a.u. Bottom: Enhanced 
and non-enhanced magnetic moment on the Rh-site in the FM 
(uDLM) state, as a function of average normalized magnetic 
moment on Fe-sites, rhpe. 


moment which are related to the magnetic disorder in 
the system. The magnetic moment of Rh is induced by 
a spin dependent hybridization of its electronic states 
with the electronic states of surrounding Fe atoms. This 
hybridization plays a crucial role during the transition 
from the FM to the AFM state. To demonstrate the 
strong covalent character of the Rh magnetism, SCF cal¬ 
culations have been performed by suppressing for the 
spin-dependent part of the exchange-correlation poten¬ 
tial (S“ = 0) that is responsible for an enhancement of 
the spin magnetic moment induced by the hybridisation 
with the Fe states. Figure [3] (open symbols) shows that 
the non-enhanced Rh magnetic moment is only ^ 25 % 
smaller than the proper one. The same is observed for 
the total energy. This demonstrates the significant role 
of the spin dependent hybridization for the formation of 
a large magnetic moment on the Rh site. As a result, 
the varying magnetic disorder in the Fe sublattice in the 
presence of the weak Rh exchange enhancement leads to 
an almost linear change of Mri, as a function of mpe. 

To investigate the stability of the FM and AFM or¬ 
dered magnetic states at finite temperature the exchange 
coupling parameters Jij have been calculated for different 
reference states: AFM, FM and DLM. These interactions 
can be seen to map the magnetic energy of the system 
onto the Heisenberg Hamiltonian accounting for the bi¬ 
linear interatomic exchange terms. The corresponding 
results are presented in Fig. IH 

Discussing these results, it is convenient to distinguish 
between the two Fe sublattices with opposite directions 
of the magnetic moments in the AFM state. For each 
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Fe atom its first and third neighbor in the Fe subsystem 
belongs to another sublattice. One can see that for all 
reference states the exchange couplings with these neigh¬ 
bours are negative indicating the trend towards the for¬ 
mation of AFM order. The interaction with the third 
Fe neighbor depends only weakly on the reference state, 
while the interactions with the first neighbour are close 
to 0 meV in the case of FM reference state and is about 
—8.0 meV for the AFM state. Since the Rh magnetic 
moment in the AFM state is equal to 0 fiB, it does not 
contribute to the magnetic energy. As a consequence 
only the Fe-Fe exchange interactions are responsible for 
the stabilization of this state. In contrast to this situ¬ 
ation, the FM order in the system can be stabilized by 
a rather strong Fe-Rh interaction since the Rh magnetic 
moments are non-zero, giving a negative contribution to 
the magnetic energy competing with the positive one due 
to the Fe-Fe interatomic exchange. Thus, the transition 
from the AFM to FM state is essentially a result of the 
competition of these interactions. 

Thus, the behaviour of the FM and AFM energy varia¬ 
tion shown in Fig. [3] and demonstrating the stabilization 
of the FM state upon increase of the magnetic disorder, 
can be attributed to the modification of the Fe-Fe ex¬ 
change coupling parameters, in particular, to a strong 
decrease of the AFM interactions with the first neigh¬ 
bors. In the DLM state the situation is very different. 
The Fe-Fe exchange interactions are rather close to those 
obtained for the FM state, but the Rh magnetic moments 
are equal to zero and therefore give no contribution to the 
magnetic energy. 

Two remarks concerning the approximations used for 
the calculations are in due place concerning the Monte 
Carlo simulations. The conventional Heisenberg Hamil¬ 
tonian should be generalized beyond the classical form: 
(i) in order to account for the contribution due to in¬ 
duced Rh magnetic moments; (ii) to account for not only 
bi-linear terms of the magnetic energy expansion but also 
for terms of higher order. The second generalization is 
required since the insufficient conventional dipolar form 
of the Hamiltonian gets appreciable corrections from bi¬ 
linear exchange coupling parameters calculated for differ¬ 
ent reference states. The general form of the expansion 
of magnetic energy around a reference state may therefor 
be represented as follows (see, e.g. [34-36|l 


E = E,^, + AEim)-J24J\e^■e,) 


( 1 ) 


i/=2 ij 


ijkl 


(ej ■ ek)iei ■ e^) -b (e; • ei){ej ■ efc)] - ... 


where AE{\Mi\) is the change in energy due to the change 



FIG. 4. (a) The Fe-Fe interatomic exchange coupling param¬ 
eters calculated for the FM (open circles), AFM (open dia¬ 
monds) and DLM (squares) states; (b) the Fe-Rh exchange 
coupling parameters calculated for the FM state (the Rh mag¬ 
netic moment is m « I^b) 


of absolute values of local spin magnetic moments. This 
can be reduced to the conventional form of the Hamilto¬ 
nian with redefined bilinear exchange interaction param¬ 
eter J 




E 

Pe:i,j 


Rh:k 


(M,.M,)(2) 


The first term in Eq. m characterizes the Fe-Fe 
transverse-fluctuation exchange energy with z, j indicat¬ 
ing sites on the Fe-sublattice and in the generalized 
form become dependent on the average magnetic mo¬ 
ment in the system. The second term describes the en¬ 
ergy changes related to longitudinal spin fluctuations on 
the Rh atoms [l3, 371 with k numbering Rh sublattice 
sites. The magnetic moments on the Fe site are denoted 
as M i(j). Considering the DLM state as reference state 
the dependence on the average magnetic moment in the 
first term can be considered in linear approximation to 
have the following form 




tFM/AFM _ tDLM 

^ij ^ij 


( 3 ) 


The term characterizing the longitudinal contribution 
was discussed previously (3^ . 

The response function Xkj occuring in Eq. ([3]) describes 
the Rh magnetic moment induced by surrounding Fe 
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atoms and is dependent on the orientation of their mag¬ 
netic moments. A linear approximation expressed by 
a constant Xkj was used, that is based on the results 
above showing the almost linear dependence of induced 
Rh magnetic moment on the average magnetic moment in 
Fe subsystem (see Fig. [3]). As the Rh magnetic moment 
occurs essentially due to the spin-dependent hybridiza¬ 
tion of the Rh electronic states with the electronic states 
of neighboring Fe atoms, it is represented in MC simula¬ 
tions through the average magnetic moment on the first 
Fe neighbor shell around the Rh atoms, leadin g to an 
approximate form for the susceptibility function 371 
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Em,. 


( 4 ) 


where the summation is performed over the magnetic mo¬ 
ments Mj corresponding to Fe atoms within the first- 
neighbor shell around the ’non-magnetic’ Rh atom on 
site i = 0. 

It should be mentioned in addition that the same DLM 
reference state for both the re-scaled FM and AFM ex¬ 
change interactions has been used. This means that 
the exchange interactions should change abruptly at the 
metamagnetic transition point that accounts for latent 
heat connected with the first-order phase transition. 

To investigate the finite temperature magnetic proper¬ 
ties of FeRh-based systems, the Monte Carlo (MC) sim¬ 
ulations have been performed based of the model Heisen¬ 
berg Hamiltonian, Eq. ([3]). In a first step the calculations 
have been performed for the pure FeRh compound. The 
magnitude of the Fe magnetic moments have been fixed 
and only changes in orientation have been considered. 
On the other hand, the magnetic moments treated as in¬ 
duced magnetic moment according to Eq. o, change 
their absolute value as well as the orientation depending 
on the orientations of the magnetic moments of the sur¬ 
rounding Fe atoms. At the same time, the total magnetic 
moment in the system can be rather small approaching 
0 at low (AFM state) and high (PM state) temper¬ 
atures. Figure [S] shows the relative magnetization as a 
function of the temperature in comparison to experimen¬ 
tal results. The calculated AFM-FM transition occurs 
at T = 320 K, rather close to the experimental value 
T = 350 K. As it was discussed above, it is caused by the 
increasing magnetic disorder in the system when the tem¬ 
perature increases. Two mechanisms are the major driv¬ 
ing force for the transition. Firstly, the disorder-induced 
modification of the exchange coupling parameters. Sec¬ 
ondly, the increase of the amplitude of randomly oriented 
fluctuations of the Rh magnetic moments in the AFM 
state due to increasing temperature-induced short-range 
FM order in the Fe subsystem (see Fig. EKc)). There¬ 
fore, it is the occurrence of magnetic moments on the Rh 
sites above a certain temperature that leads to a stabi¬ 
lization of the FM order in the system (see Fig. [D(d)). 
A further temperature increase results in a decrease of 




(d) 


FIG. 5. Temperature dependent relative magnetization 
M{T)/Mo {Mo is the magnetization at T = 0 K) obtained 
within the MC simulations: (a) M{T)/Mo for pure FeRh 
(circles) in comparison with the experimental result; squares 
and diamonds represents the results for FeRh with 5 (dia¬ 
monds) and 10 % (squares) substitution of Fe by Ni atoms; 
(b) M{T)/Mo calculated for FeRh with 1 and 2 % of Fe (closed 
symbols) in the Rh sublattice, in comparison with the results 
for FeRh (open symbols), (c) and (d) show the induced mag¬ 
netic moments on the Rh sublattice of pure FeRh at T = 200 
and 340 K, respectively. 


the Rh magnetic moment, and to a transition to the PM 
state at T = 720 K. One has to stress the asymmetry of 
the metamagnetic transition in FeRh upon heating and 
cooling of the sample. This was demonstrated recently 
by a robust experimental investigation on the formation 
of FM and AFM phases. The authors concluded that 
the formation of the AFM phase upon a temperature de¬ 
crease is dominated by a nucleation at defects in contrast 
to the formation of the FM phase for increasing temper¬ 
ature due to heterogeneous nucleation at different sites 
[ 3 ^ . The latter results are in line with the present results 
of MC simulations (e.g. see Fig. [nKc),(d)). However, dif¬ 
ferent mechanism of nucleation upon cooling requires fur- 
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ther generalization of the Hamiltonian and more sophys- 
ticated spin dynamics simulations to reproduce the tem¬ 
perature dependent behaviour of magnetization in this 
case. 

To investigate the influence of impurities on the meta- 
magnetic transition, the calculations have been per¬ 
formed for the FeRh systems with 5 at.% and 10 at.% 
substitution of Fe by Ni atoms, and with 1 at.% and 2 
at.% substitution of Rh by Fe atoms. 

The presence of Ni impurities results in a decrease of 
the temperature of the metamagnetic transition. As one 
can see in Fig. [5] (a), 5 and 10 atomic percent of Ni 
in the Fe sublattice leads to transition temperatures of 
Tm = 230 and 180 K, respectively. The decrease of Tm is 
mainly governed by the difference in the Fe-Ni exchange 
interactions when compared to the Fe-Fe exchange in¬ 
teractions. As one can see in Fig. [6l the Fe exchange 
interaction with the Ni atom at the first-neighbor posi¬ 
tion becomes positive. The Fe-Ni exchange interactions, 
when Ni occupies third-neighbor position, are negative 
but are much smaller in magnitude when compared to 
the Fe-Fe interactions. Both of these effects lead to a 
stabilization of the FM state, and as a consequence to 
a decrease the temperature of metamagnetic AFM-FM 
transition. 



FIG. 6. (a) The Fe-Ni interatomic exchange coupling param¬ 
eters calculated for the FM (open circles) and AFM (open 
diamonds) FeRh compound with 5 % substitution of Fe by 
Ni atoms: The Fe and Ni atoms correspond to different Fe 
sublattices having in the AFM state an opposite orientation 
of the magnetization. 

Substitution of 1 and 2 atomic % of Rh by Fe atoms 
results in a decrease of the transition temperature from 
Tm = 320 to Tm = 260 and 220 K, respectively. In 
contrast to Fe substitution by Ni, the decrease of Tm is 
controlled by strong FM interactions between Fe atoms 
in the different (Fe and Rh) sublattices (see Fig. [7]). 

Thus, in line with experiment, for both types of impu¬ 
rities we have obtained a decrease of the temperature of 
the metamagnetic transition. On the other hand, the ef¬ 
fect of impurities is much weaker than observed in exper¬ 
iment. This is clearly the result of approximations used 
in our calculations, in particular, for the exchange cou¬ 


pling parameters: (i) we use here the re-scaled bi-linear 
exchange interactions in the model Hamiltonian, Eq. 

(ii) the first-principles calculations of are performed 
for the collinear magnetic state at T = 0 K, which can be 
crucial for such a delicate system as FeRh. This problem 
can be avoided for example by the self-consistent DLM 
approach by Staunton et al. [l2|, that leads, however, to 
much more time-consuming calculations of the tempera¬ 
ture dependent properties. 



FIG. 7. (a) The Fe-Fe interatomic exchange coupling param¬ 
eters calculated for the FM (open circles) and AFM (open 
diamonds) FeRh compound with 1 % substitution of Rh by 
Fe atoms: one Fe atom belong to the Fe sublattice, and an¬ 
other one to the Rh sublattice. 


SUMMARY 

To summarize, we have studied here the AFM-FM 
metamagnetic transition in FeRh on the basis of the first- 
principles DFT calculations. The temperature dependent 
stability of these phases was investigated performing to¬ 
tal energy calculations for the systems with a different 
degree of magnetic disorder treated within the uDLM 
approach. The first-principles calculations supply in ad¬ 
dition the parameters (element projected magnetic mo¬ 
ments, exchange coupling parameters) for the extended 
Heisenberg model Eq. (ED. Based on this Hamiltonian, 
Monte Carlo simulations have been performed. The re¬ 
sults of both calculations allow to identify the crucial role 
and interplay of two factors: (i) the dependence of the 
Fe-Fe exchange coupling parameters on the temperature- 
governed degree of magnetic disorder in the system; (ii) 
the Rh induced magnetic moment, also dependent on the 
magnetic disorder in the system, that stabilize the FM 
state. An important observation is the competing ef¬ 
fect of the temperature dependence of these two factors. 
Increase of disorder for rising temperature leads to a de¬ 
crease of the Rh magnetic moments and as a result to 
a decrease of Fe-(Rh)-Fe FM exchange interactions re¬ 
sponsible for lowering the energy of the FM state. On 
the other hand, the decrease of the AFM Fe-Fe 
exchange interactions (see extended Hamiltonian Eq. ED 
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together with the temperature induced magnetic disorder 
leads to a stabilization of the FM state. Suppressing the 
interplay of these two effects leads to a shift of the point 
of metamagnetic transition. This was demonstrated by 
studing the impact of impurities either on the Fe or on 
the Rh sublattices. 
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